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I. INTRODUCTION 



It is a fairly long time since the ground-state criticality of the ferromagnetic Heisenberg 
chain was studied first for the classical spin. To extend the study to the quantum spin 
system required considerable effort, though it might look rather easy. The interest was 
renewed quite recently when the universality of the system was explored for any values of 
S. 

In the nineteen-sixties, Fisher explicitly calculated the free energy, the susceptibility 
and the correlation function of the classical Heisenberg chain of length L under the open 
boundary condition. According to his results, the susceptibility diverges as at T = 
in the ferromagnetic case. The correlation function decays strictly exponentially, and the 
correlation length diverges as 1/T at T = 0. 

As for the periodic chain, Joyce obtained exact expressions for the partition function, 
the correlation function and the susceptibility of the classical Heisenberg chain. His results 
agree with those of the open chain in the thermodynamic limit L — 0, as would be expected. 

For the quantum system, on the other hand, Yamada and Takahashi employed the 
Bethe-ansatz method, and numerically obtained the free energy and the zero-field suscep- 
tibility of the 5* = 1/2 ferromagnetic Heisenberg chain in the thermodynamic limit. Their 
calculations gave critical exponents as a = —1/2, u = 1 and 7 = 2. Moreover, they found 
that the free energy and the susceptibility can be expanded in ^/T at low temperatures. 

Recently, two of the present authors numerically calculated the free energy of the 
S" = 1/2 Heisenberg chain with a magnetic field in the thermodynamic limit. We concluded 
in that work that the S" = 1/2 ferromagnetic Heisenberg chain in the thermodynamic limit 
has the same scaling function of the magnetization as that of the classical ferromagnetic 
Heisenberg chain. Moreover, we conjectured that this scaling function should be universal 
for any values of S. 

In the present paper, we find the universal finite-size scaling function for the ferromag- 
netic Heisenberg chains. Thus we conclude that the ferromagnetic Heisenberg chains with 
arbitrary 5* belong to the same universality class. 

In §2, we review the scaling function and the critical exponents of the arbitrary-S* Heisen- 
berg chain in the thermodynamic limit, in order to obtain the scaling variable. We present 
the finite-size scaling function explicitly in §3. We also show numerical confirmation of the 
universality for = 1/2 and 1. In §4, we argue the origin of the universal scaling function 
of the susceptibility. 

II. SCALING FUNCTION IN THE THERMODYNAMIC LIMIT 

In this section we discuss critical exponents and the scaling law of the magnetization of 
the arbitrary-S" ferromagnetic Heisenberg chain in the thermodynamic limit. Thus we know 
what is the scaling variable for the finite-size scaling function given in the next section. 

The spin-S* ferromagnetic Heisenberg chain is defined as follows: 
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J ^ h ^ ^ 

= --^^^i- Si+i- -^^^i^ (2-1) 

i=i i=i 

where the couphng J is positive, and h denotes a magnetic field. We consider both periodic 
and open boundary conditions. 

The magnetization of the ferromagnetic Heisenberg chain has the singularity at T = h = 
0, where the susceptibility and the correlation length diverge. At low temperatures, some 
physical quantities generally show the following scaling behaviors: 

^ ~ T'rh , (2.2) 

m ~ /i^/^ at T = 0, (2.3) 
X ~ T-\ (2.4) 
e~T-^ (2.5) 

where ^ is the correlation length and x is the linear susceptibility. 

In the following, we review the scaling relations among the above critical exponents. The 
correlation function /(r) between two spins with a length r > shows the scaling behavior 
as M 

where the angular brackets (■ ■ ■) denote the thermal average, and Tq is a characteristic 
temperature. The scaling function 0(x) has the following form: 

(2.7) 

^ for X < 1, 

with the critical exponent vj. As T decreases to zero, for r <^ i^, the correction function /(r) 
behaves as 

~ T^^-'-V-", (2.9) 



where we use the scaling behavior of ^, (p.5|) . Since the correlation function /(r) remains 
finite at T = 0, we have the following relation between the critical exponents: 

2p = uri. (2.10) 

The correlation thereby reduces to 
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/(r) ~ at T = 0. 

Furthermore, we obtain the susceptibihty from (|2.6|) as 

X = f J fir)dr 

^ ' 0(r/e)d(r/O 



(2.11) 



T 

rp2l3-U-l 



(2.12) 



Comparing eqs. (p.4|) and (p.lO|) with eq. (|2.12|) , we have the following scaling relations: 

7 = -2/? + z/ + l = z/(l-r/) + l. (2.13) 
Finally, the standard scaling argument is followed by the relations 

(2.14) 
(2.15) 



A = /? + 7, 
A 5' 



Now we show the exponent relations for ferromagnetic Heisenberg models with the T = 
singularity. First, at T = 0, the argument in Appendix yields 



S J \ S 



Si- Si 



+r 



g2 



1, 



(2.16) 



where these angular brackets (■ ■ ■) denote the expectation value at the ground state instead 
of the thermal average. Comparing eq. (|2.11|) with eq. (|2.16|) , we obtain 



r] = 0. 



(2.17) 



This result ( p.l7|) is consistent with Fisher's result [jl| for S* = oo. Then the scaling relation 
( CT ), 2/3 = U7], gives 



/3 = 0. 



(2.18) 



Next, at T = and with a magnetic field, the system should be fully magnetized; hence 
eq.(|2.3|) should reduce to 



Thus we have 



171 = 1. 



(2.19) 



(2.20) 



which is consistent with (3 = through the scaling relation (|2.15|) . Now the scaling relations 
and (g3§ give 
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A = 7 = z/ + 1. (2.21) 

The above indicates that the T = singularity of ferromagnetic Heisenberg models is de- 
scribed by the one-parameter scaling. 

We are in a position to obtain all the exponents for the one-dimensional model. We 
proposed |^ the thermodynamic-limit scaling function for arbitrary magnitude of the spin 
S as follows: 

m(T, h) = mo , (2.22) 

where 

2 44 / \ 

rhoix) = -x- —x' + O [x") . (2.23) 

The scaling function of the susceptibility is hence given by 

X(T, /.) = ^Xo (2.24) 

with 

U^) = l-^x' + 0{x'). (2.25) 
Comparing ( |2.22| ) with (p. 21), we can conclude 

A = 2 and (3 = for (2.26) 
Equations ( p^.21[ ) and ( |2.26| ) give other critical exponents for the arbitrary spin as follows: 

u = l and 7 = 2 for '^S. (2.27) 

As for the classical limit S = oo, Fisher had obtained the values (|2.27|) by solving the 
model analytically. For S = 1/2, on the other hand, Yamada and Takahashi [|^-|5[ had 
obtained the values ( p.27|) numerically. 

We comment on the critical exponent a, which denotes the singularity of the specific 
heat C in the form C ~ T". In an ordinary case, we get the following scaling law: 

a + 2/5 + 7=1. (2.28) 

In the case of S* = 1/2 of our system (|2.1|), numerical calculations [^,|[ gave a as 

a = -1/2. (2.29) 



Then the values f3 = and 7 = 2 in ( 2.26 ) and ( p^.27| ) appear to be inconsistent with { 2Z 



and ( p.29| ). This inconsistency comes from a particular term of the free energy. According 
to our numerical calculations, |^ the leading term of the free energy, T^/^, does not depend 
on the magnetic field. The magnetic properties are governed by the next term of the free 
energy, T^. 
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III. FINITE-SIZE SCALING FUNCTION 



In this section, we extend the scahng forms in the previous section to the finite-size 
system. We use the dimensionless scahng parameter L/^{= T'^L) as the scaled length. 



A. General Consideration 

The finite-size scaling form near the singularity T = h = is generally given in the 
form 



^ ' ' ^ \{T/jf' {T/jyij 



(3.1) 



We substitute the critical exponents in eq.( p.l|) with eqs.( p.26D and ( |2.27] ), and we have 

Jh J \ 



m{T, h, L) = m 



(3.2) 



T^' TL J ■ 

In the thermodynamic limit L oo, eq.( |3.2| ) becomes rfiQ given by eq.( |2.23| ): 

m(T, h, L ^ oo) fn ( ] = rfiQ ( ^ ] ■ (3.3) 



2^2 



2^2 



Since mo is conjectured to be universal, we expect that the finite-size scaling function m 
in (|3.1|) is also universal, or common to all magnitudes of S. 

The function ( p.2|) is immediately followed by the scaling function of the zero-field sus- 
ceptibility of the form 



h=0 



t^^{tl) ' 



(3.4) 



where 



x{y) = 



dm{x, y) 



dx 



(3.5) 



x=0 



The finite-size scaling function x does not depend on the magnitude of the spin S', because 
m should be independent of S. 



B. Analytical calculation in the classical case (S= oo) 

Here we consider the classical Heisenberg chain in order to obtain the finite-size scaling 
function x exactly. The classical Heisenberg model is given as follows: 
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J ^ h ^ ^ 

i=l i=l 
L L 

— >• — rij • nj_|_i — as S ^ oo, 

i=l i=l 



(3.6) 



where {ni} are vectors of length unity. 

First, we consider the periodic boundary condition. Joyce obtained the partition 
function of this model, Z^, as follows: 



ZL{K)=Y.{2l + l)\f{K), 

1=0 



(3.7) 



where K = J/T and 



TT 



(3.8) 



with Ii_^_i{K) being the modified Bessel functions of the first kind. The correlation function 
is given by B 



pen 



z 



L 1=0 



M 



L-M' 



(3.9) 



where the subscript L denotes the total number of spins. The summation of the correlation 
functions (|3.9|) gives the zero-field linear susceptibility per site, Xi the form 



9 oo 

3TxP-X^,L) = l + -— -$:(/ + !) 



MK)\f^,{K)-\i^,{K)\f{K) 
Xi+,{K) - Xi{K) 



(3.10) 



We are interested in the scaling limit of eq. (|3.10|) . Noting the scaling form ( p.4|) , we leave 
only the contributions of the scaling variable TL in the limit T — and L ^ oo. At low 
temperatures. A/ behaves as follows: 



Ao[K) — > as il — > oo, 



2K 



Vi(^) 



N 



exp 



-(/ + 1 



N' 



as K —>■ oo, 



(3.11) 
(3.12) 



where N is an arbitrary positive integer. Hence we have the susceptibility at low tempera- 
tures in the following scaling form: 



XP^"(T,L) 



2J_ 

3T2 



5^(2/ + 1) exp 



d=0 



l{l + l) TL 
2 T 



+ 01- 



(3.13) 



Comparing eq. (|3.13| ) with the scaling form eq. (|3.4|) , we obtain the finite-size scaling function 
of the susceptibility as follows: 
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2 f °° 

The asymptotic behavior of x'''^" is derived as 



/(/ + !) 
2|/ 



(3.14) 



1—3 exp 



411 



6j/ 







I as y 



as y 



0, 
00, 



(3.15) 



The hming behavior in ?/ — ^ 0, or L 00, corresponds to the first term of Xo in eq.( p.25| ). 

We show in Fig.Q the temperature dependence of the hnear susceptibihty x''*'" with 
J/ {TV) fixed. In the numerical calculation of eq.( |3.10D , we sum up the terms up to / = 100, 
which is large enough for this summation to converge. We see the clear convergence to the 
scaling limit (crosses), which are given by eq.( |3.14| ). 

Next, we study the open system. In this case. Fisher gave the correlation function 
analytically as follows: 



M 



(3.16) 



where 



u{K) = cothi^: - 



K' 



(3.17) 



Summing up the correlations, we obtain the linear susceptibility of the open system as 

1 + u 2u{l-u^) 



1 -M 



(3.18) 



Taking the scaling limit T — > and L 00 with TL fixed, we have the scaling form of the 
susceptibility in ( |3.4| ) in the form 



1 — exp — 



(3.19) 



The asymptotic behavior of this scaling function is of the forms 



-( 

J_ 

3y 



i-y) 




as y 
as y 



(3.20) 



00. 



The value in ?/ — > 0, or L ^ cxd, is the same as in (|3.15|) . We show the temperature 
dependence of the linear susceptibility ( |3.19D in Fig.^ The convergence to the scaling limit 
is clear in the open system as well. 
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C. Numerical calculations in the quantum case (S= 1/2 and 1) 



Now that we analytically obtained the finite-size scaling functions ( p.l4| ) and (|3.19|) , we 
confirm them for S = 1/2 and 1 by means of numerical calculations. In order to calculate 
the susceptibility numerically, we diagonalized the Hamiltonian exactly by the Householder 
method. We treated the systems of size up to L = 14 for S = 1/2, and up to L = 8 for 
S = 1 with both periodic and open boundary conditions. We calculated the susceptibility 
as 



X(T,L) = -^VV ^ ^ . (3.21) 




First, the numerical results for the periodic boundary condition are plotted in Figs.|^ and 
^. In these figu res, w e assumed that the leading correction to the finite-size scaling form 
is of the form JT/J, because the susceptibility for S* = 1/2 can be expanded with respect 



to \jT I J in the thermodynamic limit. (This leading correction may vanish in the 

classical limit as we observe in Fig.|l|.) We extrapolated the scaling limit quadratically, using 
the three points nearest to the ordinate. 

In Fig.|^, we summarize the extrapolated date of Figs.|^ and ^ with the scaling function 
( ^.14| ). The data for S = 1/2 and S = 1 are quite consistent with the finite-size scaling 
function for S* = oo. This fact strongly suggests that the finite-size scaling function % does 
not depend on the magnitude of the spin. 

Next, for the open boundary condition, we carried out in Figs.^and ^the same analysis 
as in Figs.|^ and The finite-size scaling function (|3.19|) and the scaling-limit data from 
Figs.^ and are shown in Fig.^ The agreement is remarkable. We conclude that the finite- 
size scaling function (|3.19|) is common to all magnitudes of the spin for the open system as 
well. 



IV. DISCUSSIONS 

We showed above that the finite-size scaling function x is universal with respect to the 
magnitude of the spin both for the periodic and open boundary conditions. The scaling 
functions are given by eqs.( |3.14D and (p.l9| ), respectively. We checked the universality for 



5* = 1 /2 and 1 numerically. 

We now naturally expect that this universality of the susceptibility originates in the 
universality of the correlation function. Then we may obtain the scaling function of the 
correlation function for arbitrary S on the basis of the classical limit. For the periodic 
boundary condition, the correlation function behaves in the scaling limit as follows: 




pen 



(m, —) ior^S. (4.1) 



By using the classical correlation function (|3.9| ), we obtain the scaling function 
/P'^'''(M, J/TL) for the periodic system as follows: 
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E£o{(^ + l)exp 


r 1(1+1)] 

I 2y J 


}[exp(-f)+exp(-^); 


E£o(2/ + l)exp 


r 1(1+1)] 

I 2j/ J 





where C,i is the correlation length of the mode /, i.e. 



{l + l)T 

For the open boundary condition, the correlation function behaves as follows: 



/°P°^(M) for'^^. 



Using eq.( p.l6|) , we have 



where the correlation length ^ is J/T. 
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Appendix.. THE CORRELATION FUNCTION AT THE GROUND STATE 



We show here the derivation of the correlation function ( p.l6| ). 

It is apparent that the Heisenberg ferromagnet has 2LS + 1 number of degenerate ground 
states: 



||5'total — LS, 5'4tal — ^ 



l = 0,l,...,2LS}. 



(A.l) 



Hereafter we express these states just as \LS — I) for brevity. One of the ground states is 
given by 



\LS)^(^\S! = S). 



i=l 



We can thereby construct the rest of the ground states as 

\LS -l)= Tf- \LS) 
for / = 0, 1, . . . , 2LS, where the lowering operators Tf are defined by 



(A.2) 



(A.3) 



{2LS-iy. 

^ n{2Lsy. 



('S'total) 



(A.4) 



with 



(A.5) 



i=l 



i=l 



The coefficient in the right-hand side of ( [A .41 ) is added for the normalization of the states 



It is easy to see that the state ( |A.2|) is an eigenstate of the correlation operator Si ■ Si+r'- 

Si ■ Si+r \LS) = \LS) . (A.6) 



Now we show that all the states ( |A.3| ) are eigenstates of the operator Si ■ Si^r as well. For 
this purpose, we first observe the following commutation relation: 



'5'total) '^1 ■ ^l+r 
_ 1 

~ 2 



+r 



Thereby we have 



T; , Si ■ Si+r 







(A.7) 



(A.8) 



for / = 0, 1, . . . , 2LS. After obtaining these commutation relations, we easily see the follow- 
ing: 
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5i ■ \LS -I) =Si- Si+rTf \LS) 
= TrSi ■ \LS) 
= S%- \LS) 

= S^\LS-l) (A.9) 

for / = 0,1,..., 2L^. 

Thus an arbitrary linear combination of the ground states (|A.1| ), 

2LS 2LS 

\ij^,) = Y.Q\LS-l) with ^|q|' = 1, (A.IO) 

1=0 1=0 

gives the expectation value 

{ipgs \Sl ■ Si+rli'gs) = S'^, (A.ll) 



or eq. ( |2.16| ) in the text. 



12 



[1] M.E. Fisher: Am.J.Phys. 32(1964)343. 
[2] G.S. Joyce: Phys.Rev. 155(1967)478. 

[3] M. Takahashi and M. Yamada: J.Phys.Soc.Jpn. 54(1985)2808. 
[4] M. Yamada and M. Takahashi: J.Phys.Soc.Jpn. 55(1986)2024. 
[5] M. Yamada: J.Phys.Soc.Jpn. 59(1990)848. 
[6] H. Nakamura and M. Takahashi: J.Phys.Soc.Jpn.63(1994)2563. 
[7] S. Takada: Prog.Theor.Phys. 63(1980)1121. 

[8] T. Sakaguchi, S. Takada and K. Kubo: J.Phys.Soc.Jpn. 54(1985)861. 
[9] S. Takada and K. Kubo: J.Phys.Soc.Jpn. 55(1986)1671. 
[10] V. Privman and M. E. Fisher: Phys.Rev. B 30(1984)322. 



13 



Figure Captions 



The temperature dependence of xT'^ j J for J /{TV) = 0.5, 1.0, 1.5, 2.0, 2.5, 3.0 and 5.0 
in the classical Heisenberg chain with the periodic boundary condition. As T/J 0, 
all the data approach to the values (crosses) which are given by the finite-size scaling 
function ( ^.14[ ) with ( |3.4| ). The solid lines are guides for the eye. 

The temperature dependence of xT'^/J ioi J/{TL) = 0.5, 1.0, 1.5, 2.0, 2.5, 3.0 and 
5.0 in the classical Heisenberg chain with the open boundary condition. As T/J 0, 
all the data approach to the values (crosses) which are given by the finite-size scaling 
function ( p.l9|) with (|3.4|) . The solid lines are guides for the eye. 



3. The ^JT/J dependence of xT^/J for J/(TL) = 0.5, 1.0, 1.5, 2.0, 2.5, 3.0 and 5.0 in 
the case of the = 1/2 ferromagnetic Heisenberg chain with the periodic boundary 
condition. We used the systems of size up to L = 14. Each solid line indicates the 
extrapolating function based on the three points nearest to the ordinate. 



4. The ^T/J dependence of x^V^ for J/{TL) = 0.5, 1.0, 1.5, 2.0, 2.5, 3.0 and 5.0 
in the case of the 5 = 1 ferromagnetic Heisenberg chain with the periodic boundary 
condition. We used the systems of size up to L = 8. Each solid line indicates the 
extrapolating function based on the three points nearest to the ordinate. 

5. The finite-size scaling function xT^/J = x^'^^^{J/TL) of the ferromagnetic Heisenberg 
chain with the periodic boundary condition, ( p.l4| ). The solid curve indicates the 
scaling function x^'^^^ analytically obtained for S = oo. The numerical data, which are 
extrapolated in Figs.^ and ||, are plotted with circles for S* = 1/2 and with crosses for 
S = l. 



The ^JT/J dependence of xT^/J for J/{TL) = 0.5, 1.0, 1.5, 2.0, 2.5, 3.0 and 5.0 in the 
case of the 5* = 1/2 ferromagnetic Heisenberg chain with the open boundary condition. 
We used the systems of size up to L = 14. Each solid line indicates the extrapolating 
function based on the three points nearest to the ordinate. 



7. The ^T/J dependence of xT^/J for J/{TL) = 0.5, 1.0, 1.5, 2.0, 2.5, 3.0 and 5.0 in the 
case of the S = 1 ferromagnetic Heisenberg chain with the open boundary condition. 
We used the systems of size up to L = 8. Each solid line indicates the extrapolating 
function based on the three points nearest to the ordinate. 

8. The finite-size scaling function x^^/ J = x°^^'^{J/TL) of the ferromagnetic Heisenberg 
chain with the open boundary condition, ( p. 191 ). The solid curve indicates the scaling 
function -^"p^^ analytically obtained for S = oo. The numerical data, which are ex- 
trapolated in Figs. I and 0, are plotted with circles for S = 1/2 and with crosses for 
S = l. 
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